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Direct analytic continuation of the self energy is used to determine the effect of antiferromagnetic 
ordering on the spectral function and optical conductivity of a Mott insulator. Comparison of several 
methods shows that the most robust estimation of the gap value is obtained by use of the real part 
of the continued self energy in the quasiparticle equation within the single-site dynamical mean field 
theory of the two dimensional square lattice Hubbard model, where for U slightly greater than the 
Mott critical value, antiferromagnetism increases the gap by about 80%. 
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I. INTRODUCTION 

Quantum Monte-Carlo (QMC) evaluations of imagi- 
nary time path integrals 1 ^! 3 -! 4 - have improved to the point 
where they constitute one of the basic techniques of 
condensed matter physics. Even for fermionic problem, 
where the sign problem precludes direct simulation, the 
theoretical developments associated with single-site^ and 
cluster^ Dynamical Mean Field Theory (DMFT) have 
enabled a powerful approximate solution in terms of a 
quantum impurity model which for not too large clusters 
is sign-free or at least has a tractable sign problem. How- 
ever, while the QMC methods have proven to be very 
powerful in the study of static expectation values, ob- 
taining dynamical information has remained challenging. 
The available techniques are based either on an Exact 
Diagonalization (ED) method)^ where the number of 
states which contribute to a given response function is so 
small that level spacings become an issue, or on the an- 
alytic continuation of imaginary-time data , 10 ! 11 ' 12 which 
involves a host of other uncertainties. 

The difficulties appear with particular force in the con- 
text of the question of whether the high-T c cuprates are 
Mott insulating materials^ Recent work has suggested 
that the value of the gap and the form of the conductivity 
in the above-gap region provide important insights into 
the physics of Mott and charge-transfer insulators. 14 ' 15 ! 16 
A question of particular interest is the change in gap 
value associated with onset of antiferromagnetic order 
in a Mott insulator. One recent paper argued in favor 
of negligible changes^ while another argued for a large 
change. 14 However, determining with precision the gap 
value in theoretical model of correlated material is not 
straightforward. Fig. Q] illustrates some of the uncertain- 
ties. It shows three estimations of the local spectral func- 
tion (many body density of states) for a theoretical model 
(described more fully below) of a Mott insulator: one ob- 
tained by an ED method and two obtained by maximum 
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FIG. 1: Paramagnetic phase spectral function of the two di- 
mensional square lattice half-filled Hubbard model with near- 
est neighbor hopping and interaction parameter U = 12t com- 
puted using single site DMFT with an ED impurity solver 
(dotted lines) or an imaginary time hybridization expansion 
continuous time QMC impurity solver followed by analytic 
continuation of the measured Green's function (solid and 
dashed lines). Due to particle-hole symmetry A(ui) = A{— u>) 
only positive frequency is shown. QMC is done at inverse tem- 
perature fit = 10, a is a parameter in analytic continuation 
procedure which will be explained below. 



entropy (MaxEnt) 1 ^ analytic continuation of imaginary 
time QMC data. While the qualitative structure of the 
three estimations appear consistent, there are significant 
differences of detail, including a factor of two in the size 
of the gap which makes it difficult to compare the theo- 
retical results to data. 

In this paper we present a critical examination of dif- 
ferent methods of determining the spectral function of 
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a Mott insulator, and apply the results to the question 
of the gap value, spectral function and optical conduc- 
tivity in the paramagnetic and antiferromagnetic phases 
of the two dimensional square lattice Hubbard model. 
We study MaxEnt analytic continuation of the Green's 
function and of the self-energy, and compare the results 
to ED calculations and to direct thermodynamic eval- 
uations of the gap. We argue that continuation of the 
self energy provides the best method of minimizing the 
broadening effect of MaxEnt procedure. The self energy 
is also needed for computation of other response func- 
tions, for example the optical conductivity. We establish 
that for U near the Mott critical value the onset of anti- 
ferromagnetism increases the gap of the half filled square 
lattice Hubbard model by about 3At relative to that of 
the paramagnetic case, thus increases the gap by about 
80%. 

The rest of this paper is organized as follows. In section 
HJJwe define the model to be studied and the methods of 
solution, discuss MaxEnt and in particular its application 
to the self-energy. Section [Till presents a detailed analysis 
of the results from analytic continuation and also sum- 
marizes the interpretation of exact diagonalization data 
in light of these results. Section [IV] discusses the optical 
conductivity and section [V] is a summary and conclusion. 

II. FORMALISM 
A. Model 

We study the two dimensional Hubbard model defined 
by the Hamiltonian 
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with e v = —2t (cosp x + cosp y ) — /i. We choose the chem- 
ical potential fi = such that the electron density n = 1 . 
We shall be interested in the imaginary part of the real- 
axis electron Green's function. 

To solve the model we employ the single-site dynamical 
mean field approximation 5 which makes the approxima- 
tion that the self energy is a function of frequency only 
£(p, w) — ► On this assumption Eq. |T]) may be 

mapped on to a quantum impurity model with param- 
eters determined by a self-consistency condition. The 
essential computational task is to solve the quantum im- 
purity model to obtain the local Green's function G\ oc 
and self energy S(w). We have used two methods: a 
recently developed hybridization expansion continuous- 
time QMC procedure^ and an ED method^. The ED 
method approximates physical response functions as a 
series of poles. The QMC methods produces estimates 
of Green's functions and self energies in imaginary time. 
One must then analytically continue the QMC results to 
obtain physically relevant real frequency quantities. To 
perform the analytic continuation we used the MaxEnt 



methods pioneered in the condensed matter physics con- 
text by Gubernatis and co-workers et al 10 ' 11 ! 12 

The qualitative behavior of the model is well under- 
stood. For n = 1 and any U > 0, the ground state is an- 
tiferromagnetically ordered, and the single particle spec- 
trum is characterized by a gap, A([/), which we would 
like to compute. For U > U C 2 ~ 12t there is a gap in 
the spectrum even if antiferromagnetism is suppressed, 
and we are also interested in the value of this gap, and 
in the behavior of the spectral functions for frequencies 
near the gap edge. 

To analyze the situation more precisely we note that 
in an insulator we expect ImG(w + i0 + ) = for \uj\ < A. 
We may write the Green's function (in general a matrix) 
as 

G(p, W ) = M-H (p)-S(p, w ))- 1 . (2) 

ImG(p, uj) ^ either when ImS(p, w) ^ or, regardless 
of the value of ImX, if the "quasiparticle equation" 

det [ul - Ho(p) - ReS(p, w)] = (3) 

is satisfied for some momentum p. 

In the single-site DMFT the paramagnetic (PM) phase 
quasiparticle equation is: 

uj — RcE(o;) = — 2t(cosp x + cosp y ). (4) 

In the antiferromagnetic (AFM) phase we have 

q-1 _ f cJ-S-|-(w) 2t(cosp x + cosp y ) \ 
y 2t (cosp x + cospy) uj — Sj.(u;) J ' 

(5) 

and the quasiparticle equation is 

(uj — Re£f (a;)) (uj — Re£j(w)) = At 2 (cosp x + cosp y ) 2 . 

(6) 

We may therefore define two gaps, Ai m s, the lowest 
frequency at which ImS ^ 0, and A qp , the lowest fre- 
quency at which the quasiparticle equation is satisfied. 
If the interaction is non-vanishing, for uj > A qp phase 
space is available for a particle to decay so that we ex- 
pect A q p > Ai m £. Empirically we found A qp < Ai m s, 
suggesting that A qp ~ Ai m s for the Hubbard model. 

We may also define a third gap A M from the depen- 
dence of the particle density n on chemical potential [i. 
This is given in terms of the momentum-integrated spec- 
tral function A a = —lmG a /i: for spin a by 

n (^) =5Z / ^/(^-M) Ax(w,Ai), (7) 

where f(ui) = 1/ (exp((3uj) + 1) is the Fermi function. If 
the spectral function changes smoothly with fj, then for 
/i only slightly larger than A we would have, at T = 

n(fi) = 1 + 2 J duA(u), n = 0) + 0(fi - A) 2 , (8) 
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so that n(n) would change from 1 when \i = A. How- 
ever, it is know n 17 ' 18 ' 19 that in the PM Mott insulat- 
ing region of the single-site DMFT, the spectral function 
changes nontrivially with chemical potential, introducing 
"in-gap" states so that n begins to differ from unity at 
fi — Apm — Aot with Aq a number of order unity. It is 
not known whether this phenomenon occurs in the AFM 
phase. The fj, at which n begins to deviate from unity 
therefore provides a lower bound on the gap in the insu- 
lating state. 



B. Analytic Continuation 



physical relevance this typically includes positivity and 
a known normalization. We use a Gaussian model func- 
tion l/(\/27rcr) exp (— x 2 /(2a 2 )) with a = 5 and we have 
checked that the results do not depend on the width of 
this Gaussian. To perform the minimization we use the 
algorithm of Ref. [12. We generate spectra using a broad 
range of a and select the spectrum corresponding to the 
most probable a, according to Ref. 1 12. by calculating the 
posterior probability of a at a given G P(a\G) which up 
to a normalization factor is: 



P(a\G) = [] 



a + Xi 



cQ(A) 



(14) 



In practice QMC generates a numerical estimate F of 
a function F(t) (F(iu> n )) defined on imaginary time (or 
Matsubara frequency). For fcrmionic correlators F is re- 
lated to a spectral function A(ui) by: 



F{t) 



duj- 



1 



A(w) 



(9) 



with — 1/T the inverse temperature. 

Analytical continuation is the inversion of Eq. ^ to 
determine A given F. Unfortunately, the matrix defined 
by exp(— tu>)/ (exp(— f3u>) + 1) is extremely poorly condi- 
tioned, with many relatively small eigenvalues which, on 
inversion, greatly amplify any errors in F (i.e. differences 
between QMC estimate F and true value F) leading to 
highly unreliable estimates of A. While various attempts 
have been made to avoid this problem, the most widely 
used one is the MaxEnt method ; 10 ' 11 ! 12 MaxEnt is based 
on defining A as the function which extremizes a cost 
functional (J [{A}] which is the sum of entropy-like (S) 
and energy-like (L) terms: 



Q[{A}] = aS[{A}]-L[{A}]. 



(10) 



In Eq. (|10p , a is a temperature- like quantity that controls 
the competition between S and L. The energy-like term 
is defined in terms of the mean square misfit between the 
proposed spectrum computed from Eq. © and the QMC 
data F, which in the matrix form is: 

L = i(F-KA) T C" 1 (F-KA). (11) 

Here C is a correlation matrix which represents the un- 
certainties (statistical and systematic) in the computa- 
tion: 



(12) 



The crucial part of the method is the entropy-like term, 
which is defined in terms of a model function m(io) as 



A(w) - m(w) - A(w) In 



A(u) 



m{uj) 



(13) 



The model function is chosen to encapsulate prior in- 
formation about the function A: in the problems of 



where A is the resulting spectrum and Xi are the eigen- 
values of {A 1 ' 2 }\J\JL\ A {A 1 / 2 }. Here {A 1 / 2 } means a 
matrix with elements ^/AiSij . 

Note that uncertainties of error bar estimation from 
binning of Monte Carlo data and neglect of off diagonal 
correlation matrix elements (explained below) may intro- 
duce errors in selected a from P{a\G). In the following 
context we shall show that a reliable determination of 
gap size does not change with reasonable variance of a. 



C. Self energy 

The feature of our work is the direct continuation of the 
self energy which is related to the full Green's function 
G and the noninteracting Green's function Go by 



£ = Go 1 



G~ 



(15) 



In the Hubbard model, the self energy has the following 
asymptotic behavior^ 



1 



1 



£«,(«<;„) = U(n^)+U 2 (n^){l-(n^)) — +0(, 

lW n (lU>n) 

(16) 

If the Hartree term U{n^ a ) is subtracted the remain- 
ing is just like the Green's function with a different nor- 
malization U 2 (n_cr)(l — (n_ CT ))- We also introduce the 
cutoff frequency oj* to avoid the error in high frequencies 
and replace the self energy with frequencies above oj* by 
its known asymptotic behavior (|16p . Thus frequencies 
above oj* will not be included in the MaxEnt procedure 
and their value should be correctly reproduced if we sat- 
isfy the normalization condition. We have also verified 
that including points at to > lu* does not change our 
result since (as discussed below) the ImS at higher Mat- 
subara frequencies have a much larger error estimation 
thus contribute much smaller weight in calculating L. 

It is important to properly treat the noise correlations 
and encode them in the correlation matrix. In our work 
we do not consider errors in Go . Thus if G has an average 
value G and a measurement error SG then by expanding 



5Yi(iuj n ) = G 2 {iuj n )5G{i(jj n ) + G 3 (iuj n )5G(iLu n y 



(17) 



4 



Therefore E has a possibly non-zero shift: 

{SE(iv n )) = G- 3 (^ n )(5G(io; n ) 2 ) + • ■ • (18) 
and a correlator 

= G- 2 {iu; n )G*- 2 (iu Jm )(SG(iLu n )SG*(iLJ m )) + --- 

Structure in G means that there are important correla- 
tions in E; in particular, because at large v, G ~ 
the high frequency fluctuations in E are large. 

We estimate G(r) from continuous time measurements 
binned into uniformly discretized bins of width At. Sup- 
pose that in imaginary time the fluctuations in G are 5- 
correlated and independent of time (this may always be 
ensured by an appropriate measurement process) 

(6G(T i )SG(T j ))=g 2 S ij . (20) 

Then (assuming t are evenly discretized on [0, [3) with 
size At) 



10 - continue Z(ico n ), a=0.6 



continue £(i(fl n ), a=0.1 
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FIG. 2: Comparison of minority spin ImE(u) computed by 
subtracting continued Gq 1 and G _1 (long and short dashed 
lines), and continuing E (thick dashed and the solid lines) for 
fit = 10, U = 12t half-filled square lattice in the antiferro- 
magnetic phase. We see that the Gq 1 — G _1 one does not 
preserve positive definiteness. 



(SG(ioj n )SG*(iu m )) =/3.g 2 (AT)^„ !Wm . (21) 

Note g 2 At is expected to be independent of bin size At. 
Then Eq. (fT8|) gives zero and Eq. (fT9]) gives 

(<5E(iw n )<SE*(iw TO )) 
^G- 2 (tto n )G*- 2 (zuj m )f3g 2 (AT)S LJn ^ m , (2 "' 

which is the correlation matrix we need. Observe that 
this means that in this case the orthogonal transfor- 
mation which diagonalizes the covariance matrix is the 
transformation to Matsubara frequencies. 



D. Continuing G(r) 

The self energy may alternatively be computed by first 
continuing Go and G, inverting the continued functions, 
and subtracting. Fig. [2] compares the imaginary part of 
self-energy computed in this way to the result obtained 
by continuing E(iw„). We present two a values for each 
computation, one chosen to be near the peak of P(a\G) 
(a = 3) and one at a somewhat larger a (a = 10). 

The quantity Gq -1 — G _1 is not guaranteed to be pos- 
itive definite and we see that an unphysical sign indeed 
occurs. The difficulties are that G and Go are small 
at high frequency so that errors in the MaxEnt proce- 
dure are amplified on inversion. Errors in the position of 
the gap edge are also amplified in E. Finally, the calcu- 
lated structures in E are too broad. Another deficiency is 
seen in the PM insulating phase (not shown here) where 
Gq" 1 — G _1 fails to correctly represent the pole at the 
chemical potential which is known to exist in E. 



III. RESULTS 
A. Method 

The DMFT calculation was performed with the hy- 
bridization expansion continuous-time QMC solver—. 
Typically more than 10 9 Monte Carlo steps are made 
in each DMFT iteration, which usually takes around one 
hour CPU time on a cluster with 40 2GHz processors. 
Special attention must be paid in the paramagnetic in- 
sulating phase: in order to resolve the pole in E(ct>) one 
need a real frequency grid which has very high resolution 
in the vicinity of the chemical potential. 

We use a values that range several orders of magnitude 
to do analytic continuation. For each given a we calcu- 
late P(a\G) at the convergence, use Kramers-Kronig re- 
lation to get ReE(w) from computed ImE(w), and do a 
momentum integral of Eq.@ to get the Green's function. 

B. PM phase 

The main panel of Fig. [3] shows the imaginary part of 
continued self-energies calculated with three a value for 
the paramagnetic insulating phase of the two dimensional 
square lattice Hubbard model. We see a clear pole near 
the chemical potential but the detailed structures in side 
bands 2t < uj < 8t vary. The inset shows P(cv|G) with 
a maximum at a = 1. For a = 10 (above the maximum 
of P(cv|G)), ImE is smooth because of the regularization 
from model function. For a — 1 (at the maximum of 
P(a\G)), more structure is observed. For a = 0.1 (below 
the maximum of P(a|G)) these detailed features are more 
pronounced. The differences in the curves give some idea 
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FIG. 3: Main panel: Imaginary part of PM self-energy for 
fit = 10, U = 12t, jj, — with several a's. Inset: the posterior 
probability P(a\G) as a function of a. Note the log scale of 
a. At a = 10 (above the maximum of P(a\G)) the result is 
very smooth. At a = 1 (at the maximum of P(a\G)) more 
detailed features appear, which become more pronounced for 
a = 0.1 (below the maximum of P(a\G)). 



of the uncertainties in the process. The maximum in 
posterior probability identifies the a — 1 curve as the 
preferred continuation. 

The upper panel of Fig. [4] shows the real part of con- 
tinued self energies. The a values are the same as those 
in Fig. [3l The pole at oj = in ImE implies ReE ~ l/u> 
at low frequencies. The crossing point between ReE and 
uj + At curves gives the minimal positive solution to the 
quasiparticle equation, which gives an estimated gap size 
around to = 2.1t. Turning to the lower panel which is the 
spectral function plotted at the same w-scale, we see that 
the estimate size of 2. It is consistent with all four curves: 
ED puts its first peak slightly above lo = 2.1t, and the 
non-zero structure of QMC curves below u> = 2.\t could 
be safely considered as a result of broadening in Max- 
Ent procedure. The interesting fact is that the estimate 
of gap size is robust against a reasonable variation of 
a which provides an indication that the gap estimate is 
reliable. 

Fig. [5] shows the averaged total electron density as 
a function of chemical potential at U = I2t and var- 
ious different temperatures. This does not need ana- 
lytic continuation thus provides an independent check 
of the MaxEnt results. If we assume A(ui) changes 
slowly with temperature, then thermal fluctuation gives 
(n(T)) = (n(T = 0)) + AT 2 with A a positive num- 
ber if /i is close to the lower edge of the upper Hubbard 
band. However, on the contrary, (n) increases as temper- 
ature is decreased in Fig. [5] This suggests existence of 
"in-gap" states which increase rapidly as temperature is 
reduced. Moreover, for the lowest available temperature 
fit = 20 curve we see an almost linear dependence close 
to the band edge (2t < uj < 2.5t) which is also a result of 
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FIG. 4: Upper panel: Real part of PM self-energies for 
fit = 10, U = 12t, /j, — for several values of a, along with 
quasiparticle equation ReE(oj) = ui + At. The quasiparticle 
equation has a minimum positive solution at around lo — 2. It 
Lower panel: Spectral function constructed from the contin- 
ued self energy along with the result from ED. Vertical lines 
at to — 2.1t are drawn as eye guide. We see that all curves 
are consistent with the estimated half gap size lu = 2. It. 



"in-gap" states. Extrapolation of the curve gives 

A M = l.St which, considering the presence of "in-gap" 
states, is consistent with the discussion in section iHl that 
A qp - A M = 0.3i. 



C. AFM phase 

The main panel of Fig. [6] is the continued imaginary 
part of E(a->). The inset shows the corresponding P(a\G) 
curve. We see that the a = 0.1 curve is sharper than the 
curves for a = 0.6, 1. This is very similar with what we 
found in PM case. 

The upper panel of Fig. [7] shows the minimum pos- 
itive solution to the quasiparticle equation. We get a 
half gap size around 3.8t. This is also robust against 
changing a with a variance as small as O.lt. The dash- 
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FIG. 5: Averaged total electron density as a function of chem- 
ical potential (n(/x)) for paramagnetic U — 12t square lattice 
at inverse temperatures f3t — 5, 10, 20. We see that A M = 1.8t. 
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FIG. 6: Main panel: Imaginary part of AFM minority spin 
self-energy for /3t = 10, U = 12t, fi = with several values 
of a. Inset: P(a\G) versus a curve. We again see that the 
continued self-energy has sharper features at a = 0.1 (below 
the maximum of P(a\G)) than at a = 0.6, 1 (at or above the 
maximum of P(a\G)). 



dotted line shows the Hartree-Fock mean field solution. 
The averaged magnetization produced by Hartree-Fock 
(rn) = 0.94 agrees with the QMC value, but it predicts a 
gap size 5.6i much larger than the QMC value. 

The lower panel of Fig.[7]shows the minority spin AFM 
spectral function constructed from the continued self en- 
ergy, (^.majority M = ^minority The behavior at 
the gap edge is very sharp. The sharpness comes from 
the combined effect of Fermi surface nesting and mass 
renormalization. The half gap size is consistent with 
A qp = 3.8t. As in PM phase, ED also put its first peak 
at the gap edge. At higher frequency there is a clear dif- 
ference between ED and QMC calculations: QMC has 
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FIG. 7: Upper panel: Real part of AFM self-energy for fit — 
10, U — 12t, /i = with several values of a, and the quasi- 
particle equation ReE(tj) = ui. The Hartree-Fock mean field 
prediction (HF) is also shown. The quasiparticle equation has 
a solution at uo — 3.8t. A vertical line is drawn at uo = 3.8i 
as an eye guide. Lower panel: Spectral function constructed 
from the continued self energy along with ED result. Due to 
the particle-hole symmetry A ma j ori ty(w) = A minor ity(-aj) only 
minority spin is shown. The removal peak of spectral func- 
tion has been reflected to positive frequency and shown in the 
bottom of the lower panel. The estimate of gap edge u) — 3.8t 
is consistent with the reconstructed spectral function. 



one additional peak at around 7.5t while ED has two 
peaks at 6t and 9t and in particular, the estimate of the 
upper band edge is different. As noted in section I, ED 
and QMC+MaxEnt make different approximations and 
thus produces different self energies (especially at high 
frequencies), which translates to a difference in estimate 
of upper band edge mathematically. However, at cur- 
rent stage we do not have definite statements about its 
physical origin. 

Fig. [5] shows the averaged total electron density as 
function of [i when AFM order is allowed. As in PM 
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FIG. 8: Averaged total electron density as a function of chem- 
ical potential (n(/i)) for the Hubbard model on a square lattice 
using a self consistency condition that allows for AFM order, 
at U/t = 12 and inverse temperature fit = 5, 10, and 20. 




FIG. 9: Optical conductivities constructed from analytically 
continued QMC data for half-filled Hubbard model on a 
square lattice at U = 12t and inverse temperature fit — 10. 
The PM curve has a gap at Q ~ 4.2t and the AFM curve has 
a sharp gap edge at f2 ~ 7.6t. 



phase, this provides an independent check of MaxEnt re- 
sults. For the temperature studied {fit = 10), ~ 
3.6i < A qp . This is again consistent with the discussion 
in Section [TTJ We also see that the gap will get larger 
when temperature is reduced. Within our precision we 
cannot distinguish whether "in-gap" states exist in this 
case. We see that varying the chemical potential in the 
insulating phase leads at low T to a very sharp transition 
(visible as a slope discontinuity in the curves) be- 
tween a paramagnetic metal phase and an antiferromag- 
netic insulating phase. Even at our lowest temperature 
fit = 20 the {n{pf}) curve is apparently continuous and 
is associated with a rapid (but also apparently continu- 
ous) change of staggered magnetization (not shown here) . 
Whether the transition becomes first order as T — > re- 
mains to be studied. On increasing \i from the transi- 
tion occurs at < A qp . (At fit — 10, A p ~ 3.6i and 
A qp = 3.8t.) 

IV. OPTICAL CONDUCTIVITY 

The optical conductivity can be computed using the 
Kubo formula and the minimal coupling ansatz p — > p — 
A. The dissipative part of the conductivity is then^i 

rm = — r ' — / d2p /M -/(" + ") 

a[ ' h J-oe * J (2^) 2 O (23) 

xTr [j(p)ImG(w + n,p)j(p)ImG(w,p)] , 

where the current operator is j = 6H/5p x . 

Fig. [5] compares the paramagnetic and antiferromag- 
netic phase optical conductivities. We see that PM op- 
tical conductivity has a gap ~ 4.2t and a relatively soft 
edge, while in AFM phase it has a gap of ~ 7.6t and a 
sharp edge. This is a consequence of the large change in 



gap due to antiferromagnetism. The high-T c cuprates are 
believed to be described by t « 0.38eV (this value is for 
example the average of the even and odd parity values 
quoted in the table in section 7 of Ref. l22f ). Our result 
would imply that if U in the cuprates were of the order of 
U C 2 the optical gap would be about 2.9eV, rather larger 
than the observed 2 A w 1.8eV. 



V. CONCLUSION 

To conclude, we have presented a method to find the 
gap size from QMC DMFT calculations. We first con- 
tinue the measured Y,(iu n ) to S(w), using MaxEnt with 
the correctly estimated correlation matrix, and select a 
from the peak in the posterior probability P(a\G). We 
then plot ReS(o;) and find the lowest positive solutions 
to the quasiparticle equation. Curves corresponding to 
different values of a may give slightly different estimates, 
but this variation has been found to be small. We find 
that within our numerical accuracy the gap edge is de- 
fined by the quasiparticle equation Eq. (J3J) so (at least 
within the single-site DMFT) there are no "in-gap" ex- 
cited states arising from e.g. an excitonic binding be- 
tween a particle and a spin wave. In the paramagnetic 
phase doping produces "in-gap" states; we have estab- 
lished that the shift is about 0.3i at U > U C 2- 

For the U = 12t > U C 2 Hubbard model on the square 
lattice, we found a half gap size of 2. It in the PM phase 
and 3.8t in AFM phase. Thus for U > U C 2 antiferromag- 
netic order increases the gap relative to that of the para- 
magnetic solution by about 80%. This has been qualita- 
tively noted in Ref. [23| and Ref. which applied ED and 
QMC respectively to a model with a semicircular density 
of states, but our method provides a reliable quantitative 
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result for the square lattice model. Our finding supports 
the conclusions of Ref. [LJ that U must be somewhat less 
than the Mott critical value in the cuprates. The differ- 
ence of gap size is also apparent in the calculated optical 
conductivity, in particular the AFM optical conductivity 
is remarkably sharp near the gap edge, and has a corre- 
sponding sharp feature at the upper edge of the upper 
Hubbard band. We believe this is special to the square 
lattice, arising from the perfect nesting. 
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